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Abstract 

We study the p and 4> decays into 7r + 7r~7, 7r°7r°7 and <f> into 7r°?77 
using a chiral unitary approach to deal with the final state interac- 
tion of the MM system. The final state interaction modifies only 
moderately the large momenta tail of the photon spectrum of the 
p° — > -k + tt~j decay. In the case of cf> decay the contribution to 7r + 7r~7 
and 7r°7r°7 decay proceeds via kaonic loops and gives a distribution of 
7T7T invariant masses in which the /o(980) resonance shows up with a 
very distinct peak. The spectrum found for 4> — > 7r°7r°7 decay agrees 
with the recent experimental results obtained at Novosibirsk. The 
branching ratio for <f> — > 7r°?77, dominated by the ao(980), is also in 
agreement with recent Novosibirsk results. 
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In this work we investigate the reactions p — > 7r + 7r _ 7, 7r°7r°7 and — > 
7r + 7r _ 7, 7r°7r°7, 7r r]7, treating the final state interaction of the two mesons 
with techniques of chiral unitary theory recently developed. The energies of 
the two meson system are too big in both the p and decay to be treated 
with standard chiral perturbation theory, xPT [[jj. However, a unitary cou- 
pled channels method, which makes use of the standard chiral Lagrangians 
together with an expansion of Re T~ l instead of the T matrix, has proved to 
be very efficient in describing the meson meson interactions in all channels up 
to energies around 1.2 GeV |fj, ||, |J. The method is analogous to the effective 
range expansion in Quantum Mechanics. The work of 0] establishes a direct 
connection with xPT at low energies and gives the same numerical results 
as the work of || where tadpoles and loops in the crossed channels are not 
evaluated but are reabsorbed into the L« coefficients of the second order La- 
grangian of xPT ■ A technically much simpler approach is done in 0] where, 
only for L — 0, it is shown that the effect of the second order Lagrangian 
can be suitably incorporated by means of the Bethe-Salpeter equation using 
the lowest order Lagrangian as a source of the potential and a suitable cut 
off, of the order of 1 GeV, to regularize the loops. This latter approach will 
be the one used here, where the two pions interact in s-wave. 

The diagrammatic description for the p — > 7r + 7r~7 decay is shown in Fig. 1 
In Fig. 1 the intermediate states in the loops attached to the photon, I, 
can be K + K~ or 7r + 7r~. However, the other loops involving only the meson 
meson interaction can be also K°K° or 7r°7r° in the coupled channel approach 
of§- 

For the case of 7r°7r°7 decay only the diagrams with at least one loop 
contribute, (d), (e), (/), (g), (h), ... in Fig. 1. 

The case of the decay is analogous to the p — > 7r°7r°7 decay. Indeed, 
the terms (a), (b), (c) of Fig. 1 do not contribute since we do not have direct 
— > 7tit coupling. Furthermore, there is another novelty since only K + K~ 
contributes to the loop with a photon attached. 

The procedure followed here in the cases of 7r 7r and n°r] production is 
analogous to the one used in |J. Depending on the renormalization scheme 
chosen, other diagrams can appear |J but the whole set is calculated using 
gauge invariant arguments, as done here, with the same result. The novelty 
in the present work is that the strong interaction MM —>■ M'M' is evaluated 
using the unitary chiral amplitudes instead of the lowest order used in 0. 

We shall make use of the chiral Lagrangians for vector mesons of and 
follow the lines of ref. in the treatment of the radiative rho decay. The 



2 




(8) (h) 

Figure 1: Diagrams for the decay p — > 7r + 7r~7. 

Lagrangian coupling vector mesons to pseudoscalar mesons and photons is 
given by 

£ 2 [v(i-)) = ^(V/D + ^(v^un (i) 

where is a 3 x 3 matrix of antisymmetric tensor fields representing the 
octet of vector mesons, K*, p, cj 8 . All magnitudes involved in Eq. (|l|) are 
defined in ||. The coupling Gy is deduced from the p — > 7r + 7r~ decay 
and the Fy coupling from p — > e + e~. We take the values chosen in |7j, 
Gy = 67 MeV, Fy = 153 MeV. The meson is introduced in the scheme by 
means of a singlet, u>i, going from SU(3) to U(3) through the substitution 
V^ u — > V^ u + Iz^jf--, with J 3 the 3x3 diagonal matrix. Then, assuming ideal 
mixing for the and u mesons 
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(2) 



one obtains the Lagrangian of Eq. ([I]) substituting V^„ by V^ u , given by 
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From there one can obtain the couplings corresponding to VPP (V vector 
and P pseudoscalar) and VPP^ with the Gy term or the VPP'y with the 
Fy term. 

The basic couplings needed to evaluate the diagrams of Fig. 1 are 



+ J^P (T " Gv ) P MPW^) - fc"e*( 7 )] (4) 
£ 77r+7r - = 2ep M e M (7) 

with p M , p' the 7T + ,7r _ momenta, P^, fc M the p and photon momenta and / 
the pion decay constant which we take as f n = 93 MeV. 

The vertices of Eq. (ffi) are easily generalized to the case of K + K~ . Us- 
ing the Lagrangian of Eq. ([!]), in the first two couplings one has an extra 
factor 1/2 and the last coupling is the same. The couplings for (f)K + K~ and 
c() r )K + K^ which are needed for the <fi decay are like the two first couplings of 
Eq. (f|) substituting M p by M^, e / "(p) by e M (0) and multiplying by —l/y/2. 
In addition we shall take the values Gy = 55 MeV and Fy = 165 MeV which 
are suited to the (p — > K + K~ and (f> — ► e + e~ decay widths respectively. 

The evaluation of the p width for the first three diagrams (a), (6), (c) 
of Fig. 1 is straightforward and has been done before || [10], [□]] and in @ 



following the present formalism. We rewrite the results in a convenient way 
for our purposes 
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dMi 



(5) 



where 

E£l*l 2 = ^ 2 [/ 1 + / 2 +/ 3 ] ( 6 ) 

In Eq. d^), Mi is the invariant mass of the two 7r system and 9 the angle 
between the ir + meson and the photon in the frame where the it + tt~ system 
is at rest. The quantity I\ stands for the contribution of the first diagram 
alone, Fig. 1 (a), J3 for the second and third (6), (c) and I 2 for the interference 
between the first diagram and the other two. They are given by 
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where K is the photon momentum in the p rest frame and p, k are the mo- 
menta of the meson and the photon in the rest frame of the 7r + n~ system, 
and Di, D 2 the meson propagators in the (b), (c) Bremsstrahlung diagrams, 
conveniently written in terms of Mj and 9. 

The first term of the contact term, t fryn + 7T - , in Eq. (|J) is not gauge invari- 
ant. It requires the addition of the diagrams (6) and (c) of Fig. 1 to have a 
gauge invariant set. On the other hand the second term in the contact term 
(Fy/2 — Gy part) is gauge invariant by itself. When considering final state 
interaction of the mesons this means that the Gy part of the contact term, 
diagram (d), must be complemented by diagrams (e), (/), (g) to form the 
gauge invariant set. On the other hand the Fy/2 — Gy part of the contact 
term appears in the (d) diagram which is gauge invariant by itself. 
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The technology to introduce the final state interactions is available from 
the study of <fi — > K°K 0/ -f in [12] . There it was shown that the strong t matrix 



for the M1M2 — > M[M' 2 transition factorizes with their on shell values in the 
loops with a photon attached. The same was proved for the loops of the 
Bethe-Salpeter equation in the meson meson interaction description of ||. 
On the other hand the sum of the diagrams (d), (e), (/), (g), which appears 
now with the Gy part of the contact term (diagram (a)), could be done 
using arguments of gauge invariance which led to a finite contribution for 
the sum of the loops |5|, [13], [14]]. A sketch of the procedure is given here. The 



p — > 7r + 7r 7 amplitude can be written as e )1 (p)e u ('y)T lll/ and the structure of 
the loops in Fig. 1 is such that 

T V» =ag ^ + b QHQ» + c Q^ R u + d QV K » + e K p K v (g) 

where Q, K are the p meson and photon momenta respectively. Gauge 
invariance (T^ V K V = 0) forces 6 = and d = —a/(Q ■ K). Furthermore, in 
the Coulomb gauge only the g^ v term of Eq. (|[) contributes and the coefficient 
a is calculated from the d coefficient, to which only the diagrams (e), (/), 
of Fig. 1 contribute. For dimensional reasons the loop integral contains two 
powers less in the internal variables than the pieces contributing to the g^ LU 
term from these diagrams, since the product Q V K^ is factorized out of the 
integral. This makes the d coefficient finite. Furthermore, the MM —>■ MM 
vertices appearing there have the structure as + PJ2iPl + r yJ2i m h which 
can be recast as as + (/3 + 7) J2i m f + P J2i(Pi ~ m i )• The first two terms in 
the sum give the on shell contribution and the third one the off shell part. 
This latter term kills one of the meson propagators in the loops and does 
not contribute to the d term in Eq. @. Hence, the meson meson amplitudes 
factorize outside the loop integral with their on shell values. A more detailed 



description, done for a similar problem, can be seen in fLql , following the 
steps from Eqs. (13) to (23). 



Following these steps, as done in [12, 15], it is easy to include the ef- 



fect of the final state interaction of the mesons. The sum of the diagrams 
(d), (e), (/), (g) and further iterated loops of the meson- meson interaction, 
(h), . . . , is shown to have the same structure as the contact term of (a) in 
the Coulomb gauge, which one chooses to evaluate the amplitudes. The sum 
of all terms including loops is readily accomplished by multiplying the Gy 
part of the contact term by the factor Fi(M p , Mi) 
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Fi(M p , Mi) = 1 + (5 7r + 7r -t 7r + 7r - i7r +^- + -G K + K ~t K + K - t7r + n - (9) 

where t MlM2 ^ M ^ are the strong transition matrix elements in s-wave evalu- 
ated in and Gm x m 2 is given by 



G MlM2 (Mp, Mi) = —(a - b)I{a, b) 



(10) 



/If 2 ' M 2 

with I(a, b) a function given analytically in fl2|| . The (Fy /2 — Gy) part of the 
contact term is iterated by means of diagrams (d), (h) ... in order to account 
for final state interaction. Here the loop function is the ordinary two meson 
propagator function, G, of the Bethe-Salpeter equation, T = V + VGT, for 
the meson-meson scattering and which is regularized in [[J by means of a 
cut-off in order to fit the scattering data. The sum of all these diagrams is 
readily accomplished by multiplying the (Fy/2 — Gy) part of the contact 
term by the factor 

F 2 (Mi) = 1 + GK+ir-t-K+ir- ,7T+7T- + -jGK+K-tK+K-,TT+-K~ (11) 

By using isospin Clebsch Gordan coefficients the amplitudes £miM2,m'm' can 
be written in terms of the isospin amplitudes of 0] as 



2 _ 

^TT + TT-.TT + TT- = -t^^Mj) 

tK+K-,ir+n- = -J^KR^iMi) ( 12 ) 

neglecting the small 1 = 2 amplitudes. In Eq. (|l2l), one factor \/2 for each 
7r + 7r~ state has been introduced, since the isospin amplitudes of used in 
Eq. (0) are written in a unitary normalization which includes an extra factor 
1 / v2 for each nn state. 

The invariant mass distribution in the presence of final state interaction 
is now given by Eqs. (|5|, |], [7]) by changing in Eq. (0) 
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(13) 



The p — ► 7r 7r 7 width is readily obtained by omitting the terms I 2 ,h and 
also omitting the first term (the unity) in the definition of the Fi(M p , Mj), 



F 2 (Mj) factors in Eqs. (|9D and (|ill) and dividing by a factor two the width 
to account for the identity of the particles. 

The evaluation of the </> decay is straightforward by noting that the tree 
level contributions, diagrams (a), (6), (c) are not present now, and that only 
kaonic loops attached to photons contribute in this case. Hence, the invariant 
mass distribution for — > 7r + 7r~7 is given in this case by Eq. (|5]), changing 
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(14) 

For — > 7r°7r°7 the cross section is the same divided by a factor two to 
account fot the identity of the two 7r°'s. 
For the <j) — > 7r°r/7 case we have 
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(15) 

In Fig. 2 we show dT/dK for p — > n + n 7 decay, (dT p / dK = m p dT p /MidMi) 
The dashed-dotted line shows the contribution of diagrams 1(a), (b), (c) and 
taking Fy = 0. The dashed line shows again the contribution coming from 
diagrams 1(a), (6), (c) but now considering also the Fy contributions. Fi- 
nally, the solid line includes the full set of diagrams in Fig. 1 to account for 
final state interaction and with the Fy and Gy contributions. The process is 
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Figure 2: Photon distribution, dT/dK, for the process p — > 7r + 7r _ 7 as a 
function of the photon momentum. Solid line: spectrum including final state 
interaction of the two mesons and the Fy and Gy contributions; dashed line: 
spectrum including only the tree level diagrams (a), (5), (c) of Fig. 1 and the 
Fy and Gy contributions; dashed-dotted line: spectrum including only the 
tree level diagrams (a), (b), (c) of Fig. 1 and taking Fy = 0. The experimental 
data taken from [16| are normalized to our results. 



infrared divergent and we plot the distribution for photons with energy big- 
ger than 50 MeV, where the experimental measurements exist |16[ . We have 
also added the experimental data, given in [16| with arbitrary normalization, 
normalized to our results. 

As one can see in Fig. 2, the shape of the distribution of photon momenta 
is well reproduced. For the total contribution we obtain a branching ratio to 
the total width of the p 

B(p° -> 7T + 7r"7) = 1.18 10~ 2 for K > 50 MeV (16) 

which compares favourably with the experimental number [16|, B cxp (p° — > 
7^-7) = (0.99 ± 0.04 ± 0.15) HT 2 for K > 50 MeV. 

The changes induced by the Fy term found here reconfirm the findings of 
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Figure 3: Distribution dB/dMj for the decay — ► 7T 7r 7, with Mj the 
invariant mass of the 7r 7r system. Solid line: our prediction, with FyGy > 0. 
Dashed line: result taking FyGy < 0. The data points are from [|17] and 
only statistical errors are shown. The systematic errors are similar to the 
statistical ones jl7). The distribution for — > 7r + 7r~7 is twice the results 
plotted there. 



0. The effect of the final state interaction is small and mostly visible at high 
photon energies, where it increases dT/dK by about 25%. The branching 
ratio for B(p° — > 7r°7r°7) that we obtain is 1.4 10~ 5 which can be interpreted 
in our case as p° — ► 7o~(7r 7r°) since the tc°tt interaction is dominated by the 
a pole in the energy regime where it appears here. This result is very similar 
to the one obtained in |J. In the case one considers FyGy < 0, the result 
obtained is 1.0 10~ 4 . The measurement of this quantity may serve as a test 
for the sign of the FyGy product. 

As for the <fi — > mc-f decay, as we pointed above the <fi — > 7r + 7r~7 rate 
is twice the one of the — > 7r°7r°7. We have evaluated the invariant mass 
distribution for these decay channels and in Fig. 3 we plot the distribution 
dB j dMj for — > 7r°7r°7 which allows us to see the — > /07 contribution since 
the fo is the important scalar resonance appearing in the K + K~ — ► 7r + 7r~ 
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amplitude 0. The solid curve shows our prediction, with FyGy > 0, the 
sign predicted by vector meson dominance || . The dashed curve is obtained 
considering FyGy < 0. We compare our results with the recent ones of the 
Novosibirsk experiment [[17). We can see that the shape of the spectrum is 
relatively well reproduced considering statistical and systematic errors (the 
latter ones not shown in the figure). The results considering FyGy < are 
in complete disagreement with the data. 

The finite total branching ratio which we find for — > 7r + 7r~7 is 1.6 1CT 4 
and correspondingly 0.8 10~ 4 for the — > 7r°7r°7. This latter number is 
slightly smaller than the result given in [|I7j], (1.14±0.10±0.12) 10~ 4 , where 
the first error is statistical and the second one systematic. The result given 
rj] is (1.08 ± 0.17 ± 0.09) 10~ 4 , compatible with our prediction. The 



m 



branching ratio measured in |2(J for — > 7i + n 7 is (0.41 ±0.12 ±0.04) 10 4 . 

The branching ratio obtained for the case — > 7r°?77 is 0.87 10 -4 . The re- 
sults obtained at Novosibirsk are |9| (0.83±0.23) 10~ 4 and |L8| (0.90±0.24± 
0.10) 10 -4 . The spectrum, not shown, is dominated by the a® contribution. 

The contribution of — * / (7r + 7r~)7, obtained by integrating dT^/dMj 
assuming an approximate Breit-Wigner form to the left of the fo peak, gives 
us a branching ratio 0.44 10~ 4 . As argued above, the branching ratio for 
— > 7r°7r°7 is one half of — > 7r + 7r~7, which should not be compared to 
the one given in [17] since there the assumption that all the strength of the 
spectrum is due to the f excitation is done. As one can see in Fig. 3, we 
find also an appreciable strength for — > (77. 

We should also warn not to compare our predicted rate for — > 7r + 7r~7 
directly with experiment. Indeed, the experiment is done using the reaction 
e + e~ — > — > 7T + 7r _ 7, which interferes with the p contribution e + e~ — >• p —>■ 



7T + 7r~7 at the tail of the p mass distribution in the mass region ||21|| . Also 
the results in ^(J are based on model dependent assumptions. For these 
reasons, as quoted in |Tj|, the 7r°7r°7 mode is more efficient to study the nn 
mass spectrum. 

Our result for — > 7r°7r°7 is 50 % larger than the one obtained in owed 
to the use of the unitary K + K~ — > 7r°7r° amplitude instead of the lowest 
order chiral one. The shape of the distribution found here is, however, rather 
different than the one obtained in [Q , showing the important contribution of 
the fo resonance which appears naturally in the unitary chiral approach. 

The — > /07 decay has been advocated as an important source of in- 
formation on the nature of the fo resonance and experiments have been 
conducted at Novosibirsk 2^] and are also planned at Frascati [23], trying 
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to magnify the signal for / production through interference with initial and 



final state radiation in the e + e — > — > /o( 7r+7r )7 reaction [21], |24], |25 



The completion of the experiments [18], [H], is a significant step for- 
ward. 

Present evaluations for — > /o7 ~~ * ar e based on models assuming a 
molecule for the fo [26] with a branching ratio 1-2 10~ 5 , a gg structure 



with a value 5 10~ 5 and a gggg structure with a value 2.4 10~ 4 

The model for — > /o7 assumed in Fig. 1 is similar to the one of p7| where 
the production also proceeds via the kaonic loops. There a KK molecule is 
assumed for the fo resonance while here the realistic KK — > 7T7T amplitude 
of 0| is used. Emphasis is made in the importance of going beyond the 
zero width approximation for the resonance in P^, ffijfl . Our approach auto- 
matically takes this into account since the KK — > titt amplitude correctly 
incorporates the width of the fo resonance 0. 

We would also like to warn that the peak of the fo seen in Fig. 3 cannot 
be trivially interpreted as a resonant contribution on top of a background, 
since there are important interference effects between the fo production and 
the a background. The strength of the peak comes in our case in about 
equal amounts from the real and the imaginary parts of the amplitude for 
the process. 

The agreement found between our results for the — > 7r°7r°7, — > 7r°?77 
and experiment provides an important endorsement for the chiral unitary 
approach used here. Improvements in the future, reducing the experimental 
errors, should put further constraints on avalaible theoretical approaches for 
this reaction. 
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